Localization of phonons in ion traps with controlled quantum disorder 
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We show that the vibrations of a chain of trapped ions offer an interesting route to explore the physics of 
disordered quantum systems. By preparing the internal state of the ions in a quantum superposition, we show 
how the local vibrational energy becomes a stochastic variable, being its statistical properties inherited from the 
underlying quantum parallelism of the internal state. We describe a minimally-perturbing measurement of the 
resonance fluorescence, which allows us to study effects like Anderson localization without the need of ground- 
state cooling or individual addressing, and thus paves the way towards high-temperature ion experiments. 
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Disorder in quantum systems leads to a variety of fascinat- 
ing phenomena. In particular, the disorder-induced localiza- 
tion of particles has been at the focus of intense theoretical and 
experimental research since its discovery by P.W. Anderson 
[1]. Experiments have been performed in a variety of setups 
involving electronic transport [2], propagation of light [3] or 
sound waves [4], and disordered Bose-Einstein condensates 
[5, 6]. So far, any experiment has limitations in the degree of 
control of particle interactions or the statistical properties of 
disorder. Thus, it remains elusive to find realizations of most 
challenging situations, such as Anderson localization in the 
presence of correlated disorder or controlled interactions. In 
this work we show that the vibrations of an ion crystal can 
be used to study the physics of Anderson localization. The 
ability to control the interactions and measure the state of 
trapped ions [7], makes this setup an ideal quantum simula- 
tor of many-body systems (see e.g. [8-11]), as recently shown 
in experiments [12, 13]. Our proposal relies on the description 
of the radial vibrations of an ion chain in terms of weakly cou- 
pled harmonic oscillators [9] . These vibrations are coupled by 
means of lasers to the internal states of the ions, represented 
by a set of two-level systems, or effective spins. We employ 
a state-dependent Stark-shift that shifts the local ion trapping 
potential, to couple spins and phonons. This composite sys- 
tem can be described effectively as two different species in- 
teracting in a chain. Since the spin species is the slowest, it 
behaves like a frozen background for phonons, and makes the 
local phonon energy a true stochastic variable whose statisti- 
cal properties are determined by the spin quantum state. 

In particular, we show that: i) By preparing the ion internal 
state in a quantum superposition corresponding to a separa- 
ble state, disorder leads to the localization of phonons. This 
effect can be directly measured by studying the propagation 
of phonons along the chain, ii) Fundamental properties, like 
the phonon localization length, may be detected even at high 
temperatures, and without the need to resolve individual ions. 
This is achieved by measuring the resonance fluorescence. Hi) 
Our proposal can be extended to implement a variety of inter- 
esting situations, such as binary-correlated disorder, or local- 
ization in the presence of anharmonicities. 

General scheme.- We consider a chain of N ions with mass 
m and charge e trapped by harmonic potentials. The ions are 



coupled by the Coulomb interaction, and thus, their displace- 
ments around the equilibrium positions are described by col- 
lective vibrational modes. In particular, we focus on the vibra- 
tions in the radial - transverse to the chain - direction. Each 
ion also has two internal levels (|0), |1)) with energy differ- 
ence COo. Thus, the system Hamiltonian is (we set fi = 1) 

N N 

H = 0)o^aj+J^a n i^a n , Q? n = co? (1 +/3r n ) , (1) 

7=1 n=l 

where crj = |0 7 -)(0 7 -| — |ly)(l 7 -|, a\ (a n ) create (annihilate) ra- 
dial phonons in the vibrational mode n, and £l n are the radial- 
mode energies determined by the trapping frequency a\, and 
the harmonic correction Y n . The later is obtained from the di- 
agonalization of V jk = \ j - k\~ 3 (l - S jk ) - £ ¥i - \j - l\~ 3 8 jh 
such that % = Y.jk^jnVjk^kn^ where Mj n are the phonon 
wavefunctions. Finally, the ratio of the Coulomb repulsion 
to the trapping energy /3 = e 2 /mco^a 3 fulfills /3 <C 1, such 
that the ions correspond to weakly coupled oscillators [9]. We 
work under the simplification of equally spaced ions, although 
our results can be easily extended to inhomogeneous chains. 

We use a laser in the radial direction with frequency <Dl 
and wave vector k^, which is nearly resonant to a narrow 
transition between the internal states. Thus, the vibrational 
sidebands can be resolved, and the laser couples spins and 
phonons through the spatial dependence of the dipole cou- 
pling H L = ^o+jtej-ioL* +H.c, where of = |lj)(0/|. 
In the limit j8 <C 1, £0l can be tuned near resonance with 
all the modes by choosing 0)l — (Do = — &n + 8 n , and de- 
tunings 8 n <C £l n . Expressing the radial coordinate as xj = 
Hn ^jn{p>n + a D I V2m£l n , the coupling can be approximated 
by a multimode red-sideband Hamiltonian. Here, one consid- 
ers £2l <^ &n and r\ n = k^/^2mQ. n <C 1, where the carrier and 
remaining sidebands can be neglected, and the Hamiltonian in 
the interaction picture with respect to Ho, reads 

#int(0 = Y^FjnO+aneT 1 ** +H.C, F jn = l ^7] n ^ jn . (2) 

jn 

Under the condition \Fj n \ <C 8 n , we get an effective Hamilto- 
nian from the corrections to second order in \Fj n \/S n , 

#eff = #0 + X kjkn °f °k + L K J nm G j a t a tn , ( 3 ) 
jkn jnm 
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where = 2Fj n F£ n / 8 n is an effective coupling mediated by 
the vibrational modes, and K jnm = -F jn FJ m (8 n + 8 m )/28 n 8 m 
stands for the spin-phonon coupling. Eq. (3) can be recast 
in terms of local modes, aj = \Y,n^jn(^n a n + A^tfJ), with 
A t = ( it ) * ± ( % ) * • In the limit j3 < 1 , the Hamiltonian can 



s-p, 



be split into three different terms // e ff = H p +H s -\-K 
Hp = ^0) j a t j a j + Y< t jk a ) a k, 

j jk 

H s =Y< co o a j+ll J jk a ? a ki H *v = u H a j a ) a j 

j jk j 



(4) 



The last expression has a clear physical meaning: H v is 
an effective tight-binding Hamiltonian describing the har- 
monic coupling between ions in terms of phonon hopping; 
H s stands for a dipolar XY-model of spins subjected to a 
longitudinal field; and /f s _ p is the spin-phonon coupling in- 
duced by a state-dependent Stark- shift of the local ion en- 
ergy [14]. For non-resolved sidebands, this term would cor- 
respond to a spin-dependent light force that does not involve 
the phonon modes [15]. Typical values for the coupling con- 
stants in (4) depend on F = £2l?7l/2, T]l = k^/y/lmCOt, and 
5 = col - G)o + a\ (see Table I). 

Table I: Spin, phonon, and spin-phonon coupling strengths. 
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COj ~ G\ 
t = fiCQt 

u = -t 



J jk ~ | y -_*|3 
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cot ~ 10 4 kHz 
t ~ 10 2 kHz 
U ~ 10 kHz 
7- 1kHz 



The Hamiltonian in Eq.(4) corresponds to a chain with two 
species (phonons and spins), coupled by a density-density in- 
teraction. By choosing a\ ~ 10 MHz, 5 ~ 1MHz, j3 ~ 0.05, 
and F ~ 0.15, we obtain the hierarchy of coupling strengths 
t ~ 100kHz > U ~ 10kHz > J ~ 1kHz. Accordingly, the spin 
dynamics is much slower than the typical phonon timescales, 
such that spins may serve as an auxiliary system to induce a 
disordered background for phonons. Note that the separation 
in time scales is a consequence of the fact that we consider the 
limit (F / 8) 2 « 1. In this limit, induced effective spin-spin 
interactions are suppressed with respect to the phonon cou- 
plings, / = (F/8) 2 t. A similar situation has also been studied 
for atom-mixtures in optical lattices [16], whereas the effect 
of randomness in quantum networks was addressed in [17]. 

Controlled disorder.- We consider the initial state p(0) = 
Ps° ® Pp> where p s ° = | x l / s )( x l / s | is an arbitrary pure state with 
l^s) = !{,}^}IM>. The set {s} = {s h ...,s N }, with sj G 
{0, 1}, defines a given spin configuration, and C{ 5 j G C are the 
corresponding probability amplitudes. Due to the quenching 
of the spin dynamics, the phonon evolution can be approxi- 
mated by p p 0) = E w |c w |V 



^eff^p^e^eff \ where 
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H v +K 
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•■J^t jk a]a k + J^Ej({s})a]aj, (5) 

jk j 



and £j({s}) = (Oj + U(l + 2sy) is an on-site energy with two 
possible values £j = COj ±U. Due to the underlying quan- 
tum parallelism, the system explores simultaneously all real- 
izations of the on-site energy with probabilities p(£j({s})) = 
\ c {s}\ 2 - We stress that the binary randomness inherited from 
the frozen- spin background can be engineered at will prepar- 
ing different internal states |*P S ). In principle, arbitrary quan- 
tum states can be prepared by trapped-ion quantum logic [19]. 
Therefore, we obtain a versatile phonon model with diagonal 
disorder, formally equivalent to the so-called random binary 
alloy (RBA) [18], but in stark contrast, we can externally con- 
trol the disorder (see Fig. 1(a)). Note that the long-range hop- 
ping in Eq. (4) does not change qualitatively the phase dia- 
gram of the model. 

a) Random Binary Alloy (RBA) 




b) Random Dimer Model (RDM) 




Figure 1 : Disordered ion chain as a virtual two-leg ladder composed 
of a frozen-spin ensemble locally coupled to a system of hopping 
phonons. a) Uncorrelated RBA. b) Correlated RDM. 



Anderson localization in the RBA model- In the absence 
of disorder, phonon wavefunctions are extended over the 
whole ion chain. Conversely, impurities lead to interfer- 
ence phenomena of back- scattered waves that might local- 
ize the phonon wavefunctions. This effect, known as An- 
derson localization, occurs for any amount of disorder in ID 
systems [20]. In contrast to impurities that may consist of 
ions with a different mass [15], we consider here the spin en- 
semble as a source of disorder. As an illustrative example, 
we focus on the RBA model with uncorrelated disorder, and 
choose a separable initial spin state pf = with 
\+j) = (|0y) + | lj)) j \fl. The latter can be prepared by means 
of a laser or microwave field addressed to the internal transi- 
tion of each ion [7] . Note that such state leads to a flat prob- 
ability distribution of the local phonon energies in Eq. (5), 
which is characterized by p(£j({s})) = l/2 N . In order to ob- 
serve localization of phonons, we can either study the dynam- 
ics of single phonon excitations, or the spectral properties of 
a thermal distribution of phonons. Below, we describe these 
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two possible routes. 

a) Direct measurement of localization - The required ex- 
perimental techniques to create and measure the quantum state 
of motion have been demonstrated with single ions [21]. In 
this case, we also need cooling close to the vibrational ground 
state and individual addressing of ions. To create a single 
phonon in a ion on the chain, we assume an experimental 
procedure which is faster that the timescale for phonon tun- 
neling, l/t. This is in principle feasible, since sidebands can 
be resolved on time scales longer than \/(Q t . Considering a 
time scale t p for the pulse sequence required to create a sin- 
gle phonon, and typical values co t = 10MHz, and t = 100 
kHz, the required condition 1/^<^<1// may be fulfilled. 
Note that this is in agreement with time scales shown in [21]. 
The experiment should proceed as follows: i) Laser cooling 
close to the ground state, where the required mean phonon 
number n«l may be achieved by sideband cooling with- 
out individual resolution of the vibrational modes, ii) Initial- 
ization, where a single phonon localized at a given site / is 
created by a red-sideband coupling (single-ion version of Eq. 
(2)) that drives the transition |1/) |vac) p —> |0/) (g)aj|vac) p , 
where |vac) p is the phonon vacuum. Afterwards, internal 
states should be prepared in the coherent superposition |+ ; ) 
Hi) Switch on the disorder on-site phonon energies in (5), and 
wait for single-phonon propagation along the chain until a 
given final time tf. iv) Infer the phonon number by measur- 
ing the excitation probability of internal states. In the absence 
of disorder, the long-time dynamics at tf ^> max{l/£^} cor- 
responds to a quantum random walk where the phonon dif- 
fuses over the entire chain. Conversely, disorder yields a dras- 
tically different behavior as shown in the numerical results 
of Fig. 2, where phonons are confined within a localization 
length 5 « 10 sites. 
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Figure 2: Mean phonon number (nj(tf)) at tf = 10 3 //3ft>t, for a disor- 
dered chain, U = 3f/4, with N — 50 ions. A localized phonon at the 
center of the chain does not spread completely, but remains localized 
within the length-scale Inset: exponential tail of the distribution 
which gives us the localization length, ^ = A~ l w 10 <C N. 

b) Spectral measurement of localization - In the follow- 
ing we show how to measure phonon spectral properties. Our 
scheme perturbs minimally the spins, such that we remain in 
the quenched disorder limit implicit in Eq. (5). For this aim, 
we propose to measure the fluorescence from a cycling tran- 
sition, and show how to get the phonon localization length 



out of it, without neither cooling to the ground state nor in- 
dividual addressing. To describe the measurement scheme, 
we extend the theory of resonance fluorescence from a single 
atom [23] to an ion chain. A laser with wavevector k c in the 
radial direction drives a cycling transition between levels 
and \a,j), with detuning A and Rabi frequency £l c . Photons 
with momentum — k c are detected (see Fig. 3a)). We consider 
the limit r ^> A ^> Q c /2, where T is the natural linewidth of 
the cycling transition. In this limit, the excited level \aj) can 
be adiabatically eliminated and we get the following coupling 
between the ions and the fluorescence photons, 

j,q~k c 

Oj = i^e i2k ^\lj)(lj\ (6) 
= i ^£ (l+2ii7 c (a; + at)) |l;)<lj| + 0(T, c 2 ), 

where c\ (c q ) create (annihilate) photons with momen- 
tum q, g q is the dipole coupling of the cycling tran- 
sition, and 7] c = k c /\j2ma\. Eq. (6) describes pho- 
ton emission weighted by £l c /T, the probability ampli- 
tude for occupation of \aj). The fluorescence spectrum is 
given by y(co) = lim^ ^ g dt g drZji {0)(t)Oi (t + 
T ^ e i(fi)-fi%jT 9 ( see f or exam pi e [24]). Near co = g)l ± a\ the 
fluorescence is given by a sum of contributions from the dif- 
ferent vibrational modes. 

To be specific, we compute the resonance fluorescence 
spectrum from a phonon distribution at fine temperature T 
after Doppler cooling. In this limit, collective vibrational 
modes are not resolved, and thus the phonon density ma- 
trix describes a set of individual harmonic oscillators, p p = 
<2* -1 exp(— ^Y,j a ] a j)- The fluorescence spectrum is then 

,y ± (co) = stf £ |c w | 2 I^/^£W; } s(a) - col ±ni i} ), 

{s} jln 

(7) 

where = ^n, ^ := S^o(n + 1). Up to now, we have 
neglected the broadening of the spectral delta functions by the 
heating rates. This effect is discussed below. 

Let us show how Eq. (7) allows us to determine the phonon 
localization length. We focus on the center-of-mass (COM) 
mode (n = 0), which has the highest energy. In the ab- 
sence of disorder, the fluorescence COM peak intensity ful- 
fills Y,{s} \ c {s} 1 2 Lji * j*i^ jo^io = \(N+1), because JZ^ 
1 / y/N. Therefore the peak at the COM frequency scales lin- 
early with the number of ions, 5^- (co) °c NS(cq — C0l~ ^com)- 
Conversely, for any amount of disorder, wavefunctions be- 
come exponentially localized [20]. Since the wavefunction 
of phonons close to the COM frequency have now a spa- 
tial extent of the order of the localization length, §/, we ex- 
pect the COM fluorescence peak to saturate to S^-(co) °c 
%iS(cq — CQi — ^com), when N ^> £/. Thus, by studying the 
scaling of the COM peak with the number of ions in the chain, 
the saturation of the intensity will reveal the Anderson local- 
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ization length. We illustrate our method with a numerical cal- 
culation in Figs. 3(b) and 3(c), where the blue- sideband fluo- 
rescence is represented with and without disorder. Localiza- 
tion is evident already with a moderate ion number N « 20. 

We discuss now a few points on the measurement of res- 
onance fluorescence. Note that the sideband contributions to 
the resonance fluorescence spectrum have a linewidth of the 
order of the cooling/heating rates, T^ c , as shown, for exam- 
ple, in [23]. The sideband spectrum can be obtained with an 
heterodyne detection scheme, as discussed in [22] . During the 
measurement, radial motion of the ion chain may be heated or 
cooled with a rate of the order of Th/c = T\if^/T. In a typi- 
cal cycling transition we find T ^> co x . By choosing A = co x , 
the experiment is in the Doppler cooling regime, such that in 
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Figure 3: (a) Scheme of the cycling transition coupling the ion se- 
lectively \lj) to | ay). Selection rules impose the decay back to 
(b) Fluorescence S^-{co) (units of ^q - ) as a function of co (units of 
(Ot) for a chain of N = 100 and j8 = 0.05. Main figure: ordered case 
U = 0. Inset: disordered regime U/t = 0.25. (c) Scaling of the flu- 
orescence for p = 0.05. Blue circles: Ordered regime U = 0, which 
agrees with S^- (£2 com ) | (1 +Af) (red dashed line). Green squares: 
Disorder U/t = 0.25, where saturation shows phonon localization. 



the steady state the mean phonon number has typical values 
n = 10 [7]. Thus, the resonance fluorescence can be mea- 
sured with a steady phonon number within the Lamb-Dicke 
regime. Besides that, Doppler cooling leads to the broaden- 
ing of vibrational sidebands in Eq. (7). However, this effect 
does not spoil the measurement, as long as Th/c is of the order 
of the separation between vibrational modes, in which case, 
broadening will not affect significantly the scaling of photo- 
luminescence. Finally, note that decoherence of the internal 
states is allowed in our proposed experiment, since quantum 
coherence is not required for the effective spins to induce a 
disordered potential in Eq. (5). 

Localization-delocalization transition.- The localization 
paradigm is altered by the presence of strong interactions or 
correlated disorder, which lead to a variety of phases where a 
localization-delocalization transition holds. In the following, 
we show how to use our scheme to realize those situations. 

a) Correlated disorder - In addition to the uncorrected dis- 
order discussed above, entanglement in |\F S ) leads to corre- 
lated disorder in the phonon Hamiltonian. The simplest situ- 
ation consists of local vibrational energies which take on two 
possible values £y = C0j±U, with the constraint £y = £y + i for 
j odd (see Fig. 1(b)). This is known as the random dimer 
model (RDM) [25]. To induce such correlations, one should 
use quantum gates to prepare Bell-pairs |*F S ) = ® 7 =odd 
where |<D+) = (|0 y -)|0 y -+i) + |ly)|l y+ i))/x/2. Interestingly, 
this model presents a quantum phase transition between an 
insulating phase (v := t/U < 1), and a phase with extended 
phonon wavef unctions (v > 1), being the critical point at v c = 
1 for a nearest-neighbor model. Additionally, the disorder- 
averaged local density of states (LDOS) has a different be- 
havior in the two regimes [26]. Thus, such phase transition 
can be detected in the phonon LDOS pj(co) = (£„ JKj n 8{(0 - 
C0l — Qn)), by means of the fluorescence spectrum when the 
laser is focused onto a single ion S?-((ri) °< pj(co)- We em- 
phasize that ion traps not only allows an ideal realization of 
the RDM, but suggests the possibility of engineering a wide 
range of disorder-correlations that surpass real materials. 

b) Bose-glass phase. - Finally, radial anharmonicities in the 
trapping potentials can be enhanced by an additional standing- 
wave [9]. These non-linearities lead to repulsive (attractive) 
on-site interactions between the phonons whenever the ions 
sit at the maxima (minima) of the standing- wave laser. Com- 
bining these interactions with the effective Hamiltonian in 
Eq. (5), one reaches the disordered strongly correlated model 

jk j j 

where the interaction stregth U can be controlled by the laser 
parameters [9]. Accordingly, it is possible to study the fate of 
Anderson localization in the presence of disorder, where a yet 
to be observed insulating phase arises: the Bose glass [6, 27]. 
In contrast to the Mott insulator that arises due to the suppres- 
sion of density fluctuations in the strong repulsive regime, the 
Bose glass is a gapless and compressible phase, but still an 
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insulator due to disorder localization. Trapped ions are to be 
considered as an attractive platform to realize and detect these 
different phases. In fact, they are well suited to distingish ex- 
perimentally between Mott-insulating and Bose-glass phases. 
Following [28], it is possible to infer the compressibility via 
measurements of the boson distribution under modifications 
of a confining potential. For trapped ions, the phonons are 
harmonically confined to the center of the Coulomb crystal, 
which can be modified by tuning the trapping frequency or 
the ion number [9]. Under these variations, the phonon dis- 
tribution can be measured using standard techniques [7], and 
thus the compressibility inferred. 

Conclusions.- We have presented a scheme that uses ion 
crystals to measure phonon localization induced by disorder. 
We have shown that measuring the collective resonance flu- 
orescence suffices to prove Anderson localization, which can 
be surprisingly observed even at high phonon temperatures. 
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